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7.  A  FIVE  YEARS  TECHNICAL  RESEARCH  SUMMARY: 

(1)  Methods  to  Find  Primitive  Elements  in  Certain  Finite  Fields 

New  methods  are  developed  to  find  elements  of  order  2*>+^p  and  2*Sa 

2 

in  the  finite  field  GF(q  ),  where  q  =  2^-1  is  a  Mersenne  prime,  p  is  a 

prime  number,  and  n  is  a  divisor  of  2^  ^-1  [18,20,50].  Also  a  new  method 

2  2 

is  developed  to  find  primitive  elements  in  the  finite  field  GF(q  )  of  q 
elements  [30,36,50]. 

(2)  Fast  Transform  Algorithms  and  Applications 

Number-theoretic  transforms  are  developed  to  compute  one-dimensional 
cyclic  convolutions  [3-5,13-17,21,37].  One  advantage  of  such  transforms 
over  the  usual  discrete  Fourier  transform  is  it  uses  only  integer 
arithmetic.  Also  a  hybrid  algorithm  is  developed  to  perform  convolutions 
[21,  27,  40],  Moreover,  a  fast  polynomial  transform  is  developed  to 
compute  two-dimensional  cyclic  convolutions  [35,39] .  These  new  transform 
techniques  can  be  applied  to  digital  filtering  in  speech,  radar, 
communications  and  image  processing  [1,44-47,53-54]. 

(3)  Transform  Decoder  for  Correcting  Both  Errors  and  Erasures  of  the 
Reed-Solomon  Error-Correcting  Code 

A  fast  number-theoretic  transform  technique  is  developed  for  the 
encoding  and  decoding  of  Reed-Solomon  codes  [9-10,22-23,28-29,38].  This  new 
transform  technique  for  decoding  errors  and  erasures  of  a  (255,223) 

g 

Reed-Solomon  codes  over  GF(2  )  is  faster  than  the  conventional  Reed-Solomon 

decoding  technique  [31,33,38].  If  such  a  Reed-Solomon  code  is  concatenated 

with  a  Viterbi  decoded  convolutional  code  it  can  be  used  to  reduce  the 
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signal-to-noise  ratio  required  to  meet  the  specified  bit-error  rate 
for  deep  space  applications  and  military  applications  in  the  presence 
of  ECM.  This  transform  decoder  for  Reed-Solomon  codes  can  be  used 
in  multiple-user  communication  systems  [8,12,46]. 


(4)  3-D  Reconstruction  for  X-rays  and  Radar 

New  algorithms  are  developed  for  x-ray  reconstruction  [2,7,16, 

49,51] .  This  algorithm  improves  the  computational  speed  and  resolution  of 
x-ray  reconstruction.  A  generalized  |<u| -filter  is  developed  to  reduce 
x-ray  dosage  and  to  enhance  the  edge  of  reconstructured  pictures  [6,26, 
49,51] .  A  Weiner  filter  is  used  to  improve  the  image  quality  of  an  object 
in  the  3-D  reconstruction  by  a  weighted  average  of  successive  overlapping 
two-D  sections  of  the  object  [24,32,34,51].  A  computer  program  is  written 
for  a  specially  designed  head  mounted  frame  that  is  used  in  neurosurgery  [ 


Inner  Product  Computers,  Galois  Switching  Theory  and  Cryptography 


The  inner  product  computer  is  a  special-purpose  computational  unit 


intended  to  be  used  in  conjunction  with  a  general-purpose  computer  to 
perform  numerical  processing  tasks.  The  hardware  implementation  of  an 


inner  product  computer  is  designed,  analyzed  and  compared  with  a  general- 
purpose  computer  [11].  In  [19],  the  general  structure  of  the  Galois 


switching  theory  is  developed.  A  cryptographic  construct  is  developed  for 
the  use  of  computer  system  security  [43,  52] .  The  crypto-structure  can 


be  used  in  modern  multi-user  secure  communications  networks. 
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Abstract  -  A  quick  method  is  developed  to  find  an  element  or  order  2^*  p  in  the 
finite  field  GF(q2) ,  where  q  -  2P-1  is  a  Mersenne  prime.  Such  an 
element  is  needed  to  implement  complex  integer  transforms  of  length 
2^p  over  GF(q2)  where  3  k^_p+l. 
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10.  B.  Benjauthrit,  and  I.S.  Reed,  "On  the  Fundamental  Structure  of  Galois 

Switching  Functions,"  IEEE  Trans,  on  Computers,  Vol.  C-27,  No.  8,  August  1978, 
pp.  757-762. 

Abstract  -  It  is  shown  in  this  paper  that  the  fundamental  structure  of  Galois 

switching  functions  follows  naturally  from  that  of  Boolean  switching 
functions.  An  expanded  formula  for  deriving  multinomial  Galois 
switching  functions  is  provided  with  illustrations  of  its  application. 


Letters,  Vol.  14,  No.  21,  Oct.  12th,  1978,  pp.  697-698. 

• -  It 

Abstract  -  A  simple  method  is  developed  for  computing  elements  of  order  2  n, 
where  n  2P“*-1  and  2  <  k  <  p+1,  in  the  Galois  field  GF(qZ),  and 
q  =  2P-1  is  a  Mersenne  prime.  Such  primitive  elements  are  needed 
to  implement  complex  number-theoretic  transforms. 


21.  I.S.  Reed  and  T.K.  Truong,  "A  Fast  Computation  of  Complex  Convolution  Using 
a  Hybrid  Transform,"  IEEE  Trans.  Acoustics,  Speech,  and  Signal  Processing, 
Vol.  ASSP-26,  No.  6,  December  1978,  pp.  566-570. 

Abstract  -  In  this  paper  it  is  shown  that  the  cyclic  convolution  of  complex 
values  can  be  performed  by  a  hybrid  transform.  This  transform  is 
a  combination  of  a  Winograd  algorithm  and  a  fast  complex  integer 
transform  developed  previously  by  the  authors.  This  new  hybrid 
algorithm  requires  fewer  multiplications  than  any  previously  known 
algorithm. 


22.  I.S.  Reed  and  T.K.  Truong,  "A  Simple  Proof  of  the  Continued  Fraction 

Algorithm  for  Decoding  Reed-Solomon  Codes,"  Proceedings  IEE.  Vol.  125,  No.  12, 
December  1978,  pp.  1318-1320. 

Abstract  -  It  was  shown  recently  that  BCH  and  RS  codes  can  be  implemented  by 
Berlekamp's  algorithm  using  continued  fraction  approximations.  A 
simple  transparent  proof  of  Berlekamp's  algorithm  that  uses  such  a 
development  is  given  in  this  paper. 
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23.  I.S.  Reed,  T;K.  Truong,  and  B.  Benjauthrit,  "On  Decoding  of  Reed-Solomon 
Codes  Over  CF(32)  and  GF(64)  Using  the  Transform  Techniques  of  Winograd," 
National  Telecommunications  Conference,  Birmingham,  Alabama,  December  1978, 
pp.  20.4.1-20.4.7. 

Abstract  -  The  Winograd  algorithm  over  the  complex  number  field  is  modified 
to  compute  a  Fourier-like  transform  over  Galois  field  GF(2n) , 
where  n  =  5,6,  such  a  transform  can  be  used  to  encode  and  decode 
Reed-Solomon  (P.S)  codes  of  length  2n_l. 


24.  I.S,  Reed,  W.V.  Glenn,  G.M.  Chang,  T.K  .  Truong  and  C.M.  Chang,  "Wiener 

Filtering  of  Successive  Overlapping  Sections  of  an  X-Ray  Reconstruction," 

IEEE  Trans,  on  Biomedical  Engineering,  Vol.  BME-27,  No.  2,  February  1979, 
pp.  101-107. 

Abstract  -  In  this  paper,  a  technique  is  developed  to  improve  the  image  quality 
of  an  object  In  three-dimensional  reconstruction  by  a  weighted 
average  of  successive  overlapping  two-dimensional  sections  of  the 
object.  It  is  demonstrated  that  the  signal -to-noise  ratio  of  a 
two-dimensional  picture  can  be  improved  by  roughly  a  factor  of  2 
for  typical  x-ray  beam  shapes. - 


25.  I.S.  Reed,  T.K.  Truong  and  R.L.  Hiller,  "Correction  to  Fast  Mersenne  Prime 
Transforms  for  Digital  Filters,"  Proceedings  IEE,  Vol.  126,  No.  2,  February 
1979,  p.  203. 

Abstract  -  Some  of  the  results  in  Appendix  8.3  of  Reference  15  are  incorrect. 
In  this  letter,  we  give  the  correct  results. 


26.  I.S.  Reed,  W.V.  Glenn,  C.M.  Chang,  T.K.  Truong  and  Y.S.  Kwoh,  "Dose  Reduction 

in  X-Ray  Computed  Tomography  Using  a  Generalized  Filter,"  IEEE  Trans,  on  Nuclear 
Science.  Vol.  NS-26,  No.  2,  April  1979,  pp. 2904-2909. 

Abstract  -  Previously  with  the  use  of  computerized  phantoms  and  simulation, 

it  was  conjectured  that  the  x-ray  dosage  could  be  reduced  by  using 
the  appropriate  generalized  |w| -filter  for  x-ray  reconstruction. 

In  this  paper,  this  conjecture  has  been  confirmed  by  experiments 
that  use  data  from  an  EMI  5005  CT-scanner.  Explicitly  it  is  shown 
that  the  x-ray  dosage  can  be  lowered  to  1/11  of  the  standard  dosage 
while  still  maintaining  a  reconstruction  accuracy  which  is 
acceptable  for  some  applications. 


I.S.  Reed  and  T.K.  Truong,  "A  New  Hybrid  Algorithm  for  Computing  a  Fast 
Discrete  Fourier  Transform,"  IEEE  Trans,  on  Computers,  Vol.  C-28,  No.  7, 

July  1979,  pp.  487-492. 

Abstract  -  In  this  paper  for  certain  long  transform  lengths,  Winograd's 

algorithm  for  computing  the  discrete  Fourier  transform  (DFT)  is 
extended  considerably.  This  is  accomplished  by  performing  the 
cyclic  convolution,  required  by  Winograd's  method,  with  the  Mersenne 
prime  number-theoretic  transform  developed  originally  by  Rader. 

This  new  algorithm  requires  fewer  multiplications  than  either  the 
standard  fast  Fourier  transform  (FFT)  or  Winograd's  more 
conventional  algorithm.  However,  more  additions  are  required. 


I.S.  Reed,  T.K.  Truong  and  R.L.  Miller,  "Decoding  of  B.C.H.  and  R.S.  Codes 
with  Errors  and  Erasures  Using  Continued  Fractions,"  Electronic  Letters, 

16th  August  1979,  Vol.  15,  No.  17,  pp.  542-544. 

Abstract  -  Using  continued  fractions,  a  simplified  algorithm  for  decoding 

B.C.H.  and  R.S.  codes  is  developed  that  corrects  both  erasures  and 
errors  on  a  finite  field  GF(qm).  The  decoding  method  is  a 
modification  of  the  Fomey-Berlekamp  technique.  It  is  believed 
that  the  present  scheme  is  both  simpler  to  understand  and  to 
implement  than  more  conventional  algorithms. 


R.L.  Miller,  I.S.  Reed  and  T.K.  Truong,  "Simplified  Algorithm  for  Correcting 
Both  Errors  and  Erasures  of  Reed-Solomon  Codes,"  Proc.  IEE,  Vol.  126,  No.  10, 
October  1979,  pp.  961-963. 

Abstract  -  Using  a  finite-field  transform,  a  simplified  algorithm  for  decoding 
Reed-Solomon  codes  is  developed  to  correct  erasures  as  well  as 
errors  over  the  finite-field  GF(qm),  where  q  is  a  prime  and  m  is 
an  integer.  If  the  finite-field  transform  is  a  fast  transform,  this 
decoder  can  be  faster  and  simpler  than  a  decoder  that  uses  more 
conventional  methods. 


I.S.  Reed,  T.K.  Truong  and  R.L.  Miller,  "A  New  Algorithm  for  Computing 
Primitive  Elements  in  the  Field  of  Gaussian  Complex  Integers  Modulo  a  Mersenne 
Prime,"  IEEE  Trans,  on  Acoustics.  Speech  and  Signal  Processing,  Vol.  ASSP-27, 
No.  5,  Oct.  1979,  pp.  561-563. 

Abstract  -  A  new  method  is  developed^to  find  primitive  elements  in  the  Galois 
field  of  q2  elements  GF(q2),  where  q  is  a Mersenne  prime.  Such  ^ 
primitive  elements  are  needed  to  implement  transforms  over  GF(q  ) . 
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31.  R.L.  Miller,  T.K.  Truong  and  I.S.  Reed,  "Fast  Algorithm  for  Encoding  the 
(255,223)  Reed-Solomon  Code  over  GF(2S),"  Electronics  Letters,  13th  March, 
1980,  Vol.  16,  No.  6,  pp.  222-223. 

Abstract  -  A  new  scheme  for  reducing  the  numerical  complexity  of  the  standard 
Reed-Solomon  (RS)  encoding  algorithm  is  developed.  As  an  example, 
the  encoding  of  a  (255,223)  R.S.  code  (NASA  standard  for 
concatenation  with  convolutional  codes)  is  shown  to  require  75% 
fewer  multiplications  and  61%  fewer  additions  than  the  convert- ional 
method  of  computation. 


32.  I.S.  Reed,  W.V.  Glenn,  Y.S.  Kwoh,  T.K.  Truong  and  C.M.  Chang,  "X-Ray 

Reconstruction  of  the  Spinal  Cord  Using  Bone  Suppression,"  IEEE  Trans,  on 
Biomedical  Engineering,  Vol.  BME-27,  No.  6,  June  1980,  pp.  293-298. 

Abstract  -  In  this  paper,  a  new  method  is  developed  for  obtaining  an  x-ray 
reconstruction  of  the  soft  tissue  detail  of  the  spinal  canal. 

By  removing  the  dominant  effects  of  the  bony  vertebral  body 
within  the  projection  actual  clinical  data,  higher  quality  images  of 
the  residual  soft  tissues  components  can  be  reconstructed.  The 
intent  is  a  direct  visualization  of  the  spinal  cord  without  the 
need  for  water-soluble  contrast  (e.g.,  Metrizamide)  to  be 
installed  through  a  lumbar  or  cervical  puncture.  This  technique 
for  bone  suppresssion  also  has  potential  for  improving  visualization 
of  the  interior  of  the  mastoid  cavities  in  the  head. 


33.  I.S.  Reed,  R.L.  Miller  and  T.K.  Truong,  "Efficient  Program  for  Decoding  the 
(255,233)  Reed-Solomon  Code  over  GF(2°)  with  both  Errors  and  Erasures  Using 
Transform  Decoding,"  Proc.  IEE,  Vol.  127,  Pt.  E,  No.  4,  July  1980,  pp.  136-142. 

Abstract  -  To  decode  a  (255,223)  Reed-Solomon  code  over  GF(2®),  a  fast  Fourier- 
like  transform  over  GF(2°)  has  been  developed  to  compute  the 
syndromes  and  the  error-erasure  vectors  of  the  transmitted  code 
words.  This  new  simplified  transform  decoder  is  implemented  in  a 
program  on  a  digital  computer.  The  (255,223)  Reed-Solomon  code 
over  GF(2®)  is  being  proposed  as  a  NASA  standard  for  concatenation 
with  a  (7,  *s)  conventional  code.  In  a  simulation,  random  code 
words  were  corrupted  by  random  error  and  erasure  patterns,  and 
decoded  whenever  theoretically  possible.  A  matrix  of  execution 
times  for  this  new  transform  decoder  under  varying  sets  of  errors 
and  erasure  patterns  is  included  in  the  paper.  This  matrix 
demonstrates  that  the  speed  of  the  new  decoder  is  between  three  and 
seven  times  faster  than  the  software  R-S  decoder  developed  previously 
by  NASA. 


me 


34.  I.S.  Reed,  W.V.  Glenn,  Y.S.Kwoh  and  T.K.  Truong,  "A  Bandpass  Filter  for  the 
Enhancement  of  an  X-Ray  Reconstruction  of  the  Tissue  in  the  Spinal  Canal," 
IEEE  Trans,  on  Biomedical  Engineering,  Vol.  BME-27,  No.  12,  December  1980, 
pp.  736-738. 

Abstract  -  In  this  communication,  a  new  bandpass  reconstruction  filter  is 

developed  to  partially  remove  the  low  spatial  frequencies  of  the 
bone  and  the  soft  tissue  in  an  X-ray  reconstruction  of  a  lumbar 
spine.  This  partial  removal  of  the  low  frequencies  suppresses  the 
bony  vertebral  body  and  the  soft  tissue  components  within  the 
projections  of  actual  clinical  data.  It  also  has  the  effect  of 
enhancing  the  sharp  edges  of  the  fatty  tissue  surrounding  the 
spinal  cord  region.  The  intent  of  this  effort  is  to  directly 
visualize  the  spinal  cord  without  the  need  for  water-soluble 
contrast  (e.g.,  metrizamide)  to  be  installed  through  lumbar 
punctures . 


35.  T.K.  Truong,  I.S.  Reed,  R.  Lipes  and  C.  Wu,"  On  the  Application  of  a  Fast 
Polynomial  Transform  and  the  Chinese  Remainder  Theorem  to  Compute  a  Two- 
Dimensional  Convolution,"  IEEE  Trans,  on  Acoustics,  Speech  and  Signal  Processing, 
Vol.  ASSP-29 ,  No.  1,  Feb.  1981,  pp.  91-97. 


Abstract  - 


In  this  paper,  a  fast  algorithm  is  developed  to  compute  two-dimensional 
convolutions  of  an  array  of  dixd2  complex  number  points  where  d2  *  2m 
and  dj  -  2m-r+1  for  some  1  <  r  <  m.  This  new  algorithm  requires 
fewer  multiplications  and  about  the  same  number  of  additions  as  the 
conventional  FFT  method  for  computing  the  two-dimensional  convolution. 
It  also  has  the  advantage  that  the  operation  of  transposing  the 
matrix  of  data  can  be  avoided. 


36.  R.L.  Miller,  I.S.  Reed  and  T.K.  Truong,  "A  Theorem  for  Computing  Primitive 

Elements  in  the  Field  of  Complex  Integers  of  a  Characteristic  Mersenne  Prime,' 
IEEE  Trans,  on  Acoustics.  Speech,  and  Signal  Processing.  Vol.  ASSP-29,  No.  1, 
February  1981,  pp.  119-120. 


Abstract  - 


A  method  developed  previously  [30]  for  computing  primitive  elements 
in  GF(q^),  where  q  is  a  Mersenne  prime,  is  shown  not  to  generalize 
to  other  Galois  fields.  The  method  will  be  successful  in  finding 
primitive  elements  of  GF(qn)  if  and  only  if  q  is  a  Mersenne- prime 
and  n  ■  2. 


I.S.  Reed,  T.K.  Truong,  B.  Benjauthrit,  and  C.  Wu,  A  Fast  Algorithm  for 
Computing  a  Comp lex -Number  Theoretic  Transform  for  Long  Sequences,"  IEEE  Trans, 
on  Acoustics,  Speech,  and  Signal  Processing,  Vol.  ASSP-29,  No.  1,  February 
1981,  pp.  122-124. 

Abstract  -  This  correspondence  describes  an  algorithm  for  computing  a  complex- 
number  theoretic  transform  of  long  sequences.  Such  a  transform 
technique  can  be  used  to  compute  the  convolution  of  two  sequences  of 
complex  numbers.  Emphasis  is  given  to  a  transform  of  length  9*8*31 
2232.  Such  a  transform  can  be  used  to  perform  the  matched-filter 
correlation  of  raw  radar-echo  data  and  the  range  response  function 
for  producing  images  from  synthetic  aperture  radar  (SAR)  data. 


I.S.  Reed,  T.K.  Truong,  R.L.  Miller  and  J.  P.  Huang,  "Fast  Transforms  for 
Decoding  Reed-Solomon  Codes,"  IEE  Proc. ,  Vol.  128,  Pt.  F,  No.  1,  February  1981, 
pp.  9-14. 

Abstract  -  In  the  paper  it  is  shown  that  the  Chinese  remainder  theorem  when 
coupled  with  a  modification  of  Winograd's  method  can  be  used  to 
compute  Fourier-like  transforms  over  GF(2m),  where  m  *  2,3,. ..,8. 
These  new  transform  techniques  are  to  decode  Reed-Solomon  codes  of 
block  length  2m-l.  The  results  are  shown  to  be  more  efficient  than 
the  more  conventional  method. 


I.S.  Reed,  H.M.  Shao,  and  T.K.  Truong,  "Fast  Polynomial  Transform  and  its 
Implementation  by  Computer,"  IEE  Proc. ,  Vol.  128,  Pt.  E,  No.  1,  March  1981, 
pp.  50-60. 

Abstract  -  Recently  a  new  algorithm  was  developed  to  compute  two-dimensional 
cyclic  convolution  by  what  is  called  the  FPT  (fast  polynomial 
transform)  algorithm.  In  this  paper  this  new  algorithm  is  further 
studied  and  implemented  on  a  general  purpose  computer.  Methods  to 
cope  with  finite  core  memory  limitations  are  presented  and  actual 
computational  speed  performances  are  listed. 

I.S.  Reed,  T.K. Truong,  and  B.  Benjauthrit,  "Addendum  to  a  New  Hybrid  Algorithm 
for  Computing  a  Fast  Discrete  Fourier  Transform,"  IEEE  Trans,  on  Computers, 

Vol.  C-30,  No.  6,  June  1981,  pp.  453-454. 

Abstract  -  Recently,  the  authors  [27]  proposed  a  hybrid  algorithm  for  computing 
the  discrete  Fourier  transform  (DFT)  of  certain  long  transform 
lengths.  In  that  technique,  a  Wlnograd-type  algorithm  was  used  in 
conjunction  with  the  Mersenne  prime-number  theoretic  transform  to 
perform  a  DFT.  Even  though  this  technique  requires  fewer  multi¬ 
plications  than  either  the  standard  fast  Fourier  transform  (FFT) 


or  Winograd's  more  conventional  algorithm,  it  increases  the  number 
of  additions  considerably.  In  this  letter  it  is  proposed  to  use 
Winograd’s  algorithm  for  computing  the  Mersenne  prime-number 
theoretic  transform  in  the  transform  portion  of  the  hybrid  algorithm. 
It  is  shown  that  this  can  reduce  significantly  the  number  of 
additions  while  still  maintaining  about  the  same  number  of 
multiplications . 

41.  Y.S.  Kwoh,  I.S.  Reed,  T.K.  Truong  and  W.V.  Glenn,  "A  New  CT  Collimator  for 
Producing  Two  Simultaneous  Overlapping  Slices  from  One  Scan,"  IEEE  Trans .  on 
Biomedical  Engineering ,  Vol.  BME-28,  No.  9,  September  1981,  pp.  664-667. 

Abstract  -  In  this  paper,  a  modification  of  the  EMI  CT  5005  collimator  is  used 
to  simultaneously  produce  two  overlapping  two-dimensional  sections 
of  the  object  from  one  scan.  Such  a  new  collimator  not  only  achieves 
a  speed  improvement  in  the  generation  of  successive  overlapping 
sections,  but  also  yields  a  reduction  of  the  X-ray  dosage  used  by 
the  standard  EMI  CT  5005  collimator.  It  is  demonstrated  that  the 
quality  of  the  images  from  these  : two  slices  almost  equals  the 
quality  of  the  images  that  are  produced  by  the  EMI  CT  5005 
collimator  from  two  separate  scans. 


42.  H.M.  Shao,  I.S.  Reed,  T.K.  Truong  and  Y.S.  Kwoh,  "An  Improved  CT-Aided 

Stereotactic  Neurosurgery  Technique,"  Proc.  of  the  Fifth  Annual  Symposium  on 
Computer  Applications  in  Medical  Care,  November  1981,  pp.  591-595. 

Abstract  -  In  this  paper,  computer  software  is  developed  for  a  specially 
designed  head  mounted  frame  that  is  used  in  neurosurgery .  This 
system  is  applied  to  stereotactic  surgery  .  It  is  shown  to  provide 
a  simple,  fast  and  accurate  tool  for  brain  neurosurgery. 


43.  D.K.  Kravitz  and  I.S.  Reed,  "Extension  of  RSA  Crypto-Structure:  A  Galois 

Approach,"  Electronics  Letters,  18th  March  1983,  Vol.  18,  No.  6,  pp.  255-256. 

Abstract  -  The  Euler  totient  function  and  Euler-Fermat  theorem  utilized  in  the 
RSA  scheme  are  extended  from  the  integers  to  polynomials  over  finite 
fields.  The  new  scheme  is  suited  for  both  privacy  and  authentication 
implementations,  as  is  its  predecessor.  The  security  of  the  system 
rests  in  part  on  the  difficulty  of  determining  the  degrees  of  the 
irreducible  factors  of  a  high-degree  polynomial. 
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44.  T.K.  Truong,  I.S.  Reed,  C.-S.  Yeh  and  H.M.  Shao,  "A  Parallel  VLSI  Architecture 
for  A  Digital  Filter  of  Arbitrary  Length  Using  Fermat  Number  Transforms," 
Proceedings  of  1982  IEEE  International  Conference  on  Circuits  and  Computers, 
Sept.  28  -  Oct.  1,  1982,  New  York,  pp.  574-578. 

Abstract  -  In  this  paper  a  parallel  and  pipeline  architecture  is  developed  to 

compute  the  linear  convolution  of  two  sequences  of  arbitrary  lengths 
using  the  Fermat  number  transform  (FNT) .  Such  a  new  architecture 
is  developed  to  realize  the  overlap-save  method  using  one  FNT 
and  several  reverse  FNT's  of  128  points.  The  generalized  overlap- 
save  method  alleviates  the  usual  dynamic  range  limitation  in  FNT's  of 
long  transform  lengths.  Its  architecture  is  regular,  simple, 
expandable,  and  naturally  suitable  for  VLSI  implementation. 

45.  T.K.  Truong,  I.S.  Reed,  and  C.-S.  Yeh,  "A  Parallel  Architecture  for  Computing 
Cyclic  Convolutions,"  Proceedings  of  the  Trends  in  Electronics  Conference, 

Dec.  6-8,  1982,  Hong  Kong,  pp.  1-6. 

Abstract  -  In  this  paper,  a  parallel  architectural  structure  is  developed  to 
compute  one-dimensional  cyclic  convolutions.  This  structure  is 
based  on  the  Chinese  remainder  theorem  and  Kung's  systolic 
array.  To  compute  a  d-point  cyclic  convolution  the  structure 
needs  d/2  inner  product  cells.  The  utilization  of  this  structure 
is  high.  Also  it  is  simple  and  regular,  and  therefore  suitable 
for  VLSI  implementation. 

46.  T.K.  Truong,  L.J.  Deutsch,  I.S.  Reed,  I.S.  Hsu,  K.  Wang,  and  C.-S.  Yeh, 

"The  VLSI  Design  of  a  Reed-Solomon  Encoder  Using  Berlekamp’s  Bit-Serial 
Multiplier  Algorithm,"  Proceedings  of  the  Third  Caltech  Conference  on  Very 
Large  Scale  Integration.  Calif.  Institute  of  Technology,  Pasadena,  Calif., 

1983,  pp.  303-329. 

Abstract  -  E.R.  Berlekamp  has  developed  for  the  Jet  Propulsion  Laboratory 
a  bit-serial  multiplication  algorithm  for  the  encoding  of 
Reed-Solomon  (RS)  codes,  using  a  dual  basis  over  a  Galois  field. 

The  conventional  RS-encoder  for  long  codes  often  requires  look-up 
tables  to  perform  the  multiplication  of  two  field  elements. 
Berlekamp's  algorithm  requires  only  shifting  and  exclusive-OR 
operations.  It  is  shown  in  this  paper  that  the  new  dual -basis 
(255,  223)  RS-encoder  can  be  realized  readily  on  a  single  VLSI  chip 
with  NMOS  technology. 

47.  I.S.  Reed,  C.-S.  Yeh,  and  T.K.  Truong,  "A  VLSI  Architecture  for  Digital 
Filters  Using  Complex  Number-Theoretic  Transforms,"  Proceedings  1983  Interaationj 
Conference  on  Acoustics,  Speech,  and  Signal  Processing,  April  14-16,  1983, 
Boston,  Mass.,  pp.  923-926. 

Abstract  -  In  this  paper  a  parallel  architecture  is  developed  to  realize  a 
digital  filter.  First  a  systolic  array  is  used  to  compute  a 
248-point  complex  number-theoretic  transform  (CNT).  Next  an 
algorithm  is  developed  to  realize  a  digital  filter  that  uses 
248-point  CNT's  and  a  generalization  of  the  overlap-save  method. 
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This  algorithm  solves  the  conflict  between  long  transform  lengths 
and  a  wide  dynamic  range  associated  with  the  number-theoretic 
transform.  Finally  this  algorithm  is  mapped  to  a  parallel 
architecture.  This  architecture  is  simple,  regular  and  expandable, 
and,  hence,  is  suitable  for  VLSI  implementation. 

48.  I.S.  Reed,  and  T.K  .  Truong,  "New  Syndrome  Decoder  for  (n,l)  Convolutional 
Codes,"  Electronics  Letters,  Vol.  19,  No.  9,  25th  April,  1983,  pp.  344-346. 

Abstract  -  This  letter  presents  a  new  syndrome  decoding  algorithm  for  the  (n,l) 
convolutional  codes  (CC)  that  is  different  and  simpler  than  the 
previous  syndrome  decoding  algorithm  of  Schalkwijk  and  Vinck. 

The  new  technique  uses  the  general  solution  of  the  polynomial  linear 
Diophantine  equation  for  the  error  polynomial  vector  E(D).  A 
recursive,  Viterbi-like,  algorithm  is  developed  to  find  the 
minimum  weight  error  vector  t( D) .  An  example  is  given  for  the 
binary  non-systematlc  (2,1)  CC. 


B.  Dissertations  and  Abstracts 


49.  Y.S.  Kwoh,  Application  of  Finite  Field  Transforms  to  Image  Processing  and 

X-Ray  Reconstruction,  Ph.D.  Dissertation,  University  of  Southern  California, 
January  1977. 

Abstract 

This  dissertation  develops  and  utilizes  an  important  class  of  finite 
field  transforms  which  are  similar  to  the  discrete  Fourier  transform. 
Specifically  a  transform  is  defined  here  over  GF(p) ,  where  p  is  a  prime  of 
the  form  k • 2n+l ,  where  k  and  n  are  integers. 

In  this  thesis  the  above  transform  is  used  first  to  filter  a  two-dimension< 
picture.  Then  the  results  are  compared  with  the  standard  FFT. 

Another  area  of  image  processing  involves  the  reconstruction  of  an  x-ray 
image  from  one-dimensional  projection  data.  It  is  also  applied  to  3-D 
reconstruction  with  a  consequent  speed  and  accuracy  advantage  over  the 
conventional  method. 

Important  problem  of  3-D  reconstruction  which  arose  in  the  above  are 
also  studied  in  detail.  These  include  a  development  of  a  new  generalized 
|w|-filter.  The  parameters  of  the  generalized  ]u)| -filter  presented  here 
have  good  reconstruction  accuracy  as  well  as  flexibility. 

Also,  it  turns  out  that  the  development  of  the  generalized  |u>j -filter 
leads  to  a  very  impressive  result,  namely,  a  substantial  reduction  of  x-ray 
dosage.  Hence  the  combination  of  the  finite  field  transform  and  the 
generalized  |w| -filter  promise  to  yield  significant  improvements  over  existing 
systems . 
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50.  K.Y .  Liu,  Efficient  Digital  Transform  Algorithms  and  Architectures  for  Signal 
Processing  and  Error  Correcting  Codes,  Ph.D.  Dissertation,  University  of 
Southern  California,  June  1977. 

Abstract 

Several  classes  of  efficient  digital  transform  algorithms  and  architectures 
that  have  application  to  dqgltal  signal  processing  and  error-correcting  codes 
are  presented  in  this  dissertation. 

It  is  shown  that  both  high  radix  and  real-valued-input,  fast-Fourier- 
transform  (FFT)  algorithms  can  be  used  with  number- theoretic  transforms  on 
the  finite  field  GF(q^) ,  where  q  ■  2^-1  is  a  Mersenne  prime,  when  the 
transform  length  d  divides  2^+^.  Also  shown  is  a  highly  efficient  mixed-radix 
FFT  over  GF(q^)  for  d  ■  2^p  with  0<k<p+l.  This  new  FFT  algorithm  requires 
substantially  fewer  multiplications  than  the  conventional  FFT  and  the  Uinograd 
FFT  algorithm  over  the  complex  numbers. 

A  special  number-theoretic  transform  that  can  be  computed,  using  a 
high-radix  FFT,  is  defined  on  primes  of  the  form  (2n-l)2n+l.  Also  Included 
are  methods  for  finding  these  primes  and  the  primitive  dt'1  roots  of  unity 
in  the  finite  field  of  integers  modulo  such  primes. 

A  high-radix  FFT  with  generator  Y  “  3  over  GF(Fn)  where  Fn  **  22n+l  is  a 
Fermat  prime  is  found  to  be  useful  for  the  fast  encoding  and  decoding  of  Reed- 
Solomon  (RS)  codes.  A  new  class  of  highly  efficient  LSI  logic  structures, 
called  assembly-line  structures,  are  studied.  Real-time  digital  signal 
processors  can  be  constructed  by  combining  specialized  LSI  chips,  utilizing 
such  logic  structures. 


51.  C.M.  Chang.  3-D  Reconstruction  and  the  Technique  of  X-Ray  Computerized 

Tomography,  Ph.D.  dissertation,  University  of  Southern  California,  July  1979. 

Abstract 


In  this  dissertation  it  is  shown  that  the  fast  digital  convolution 
technique  used  for  parallel  x-ray  beams  can  be  used  also  to  reconstruct  a 
density  function  for  diverging  x-ray  beams.  The  convolution  method  when  used 
with  a  generalized  |m| -filter  yields  good  accuracy  and  a  flexibility  to  filter 
unwanted  noise.  It  is  shown  that  the  x-ray  dosage  can  be  lowered  to  1/11  of 
the  standard  dosage  while  still  maintaining  a  reconstruction  accuracy  that  is 
acceptable  for  many  applications. 

An  important  problem  of  x-ray  computed  tomography  is  to  improve  the  image 
quality  so  that  the  x-ray  image  can  be  diagnosed  accurately.  This  problem  is 
studied  here  in  some  detail.  A  technique  is  developed  to  improve  the  CT 
image  quality  of  an  object  in  3-D  reconstruction  by  a  weighted  average  of 
successive  overlapping  2-D  sections  of  the  object. 


Finally,  a  new  method  is  developed  to  reconstruct  the  soft  tissue  detail 
of  the  spinal  canal.  By  removing  the  dominant  effects  of  the  bony  vertebral 
body  within  the  projection  data,  it  is  demonstrated  that  higher  quality 
images  of  the  residual  soft  tissue  components  can  be  reconstructed. 

David  William  Kravitz,  The  Application  of  Finite  Polynomial  Rings  to  Number 
Theoretic  Crypto-Systems,  Ph.D.  Dissertation,  University  of  Southern  California, 
May  1982. 

Abstract 

A  cryptographic  construct  is  proposed  based  on  traditional  information- 
theoretic  and  modem  computational  complexity-theoretic  principles .  The 
methodology  developed  is  suitable  for  the  implementation  of  both  conventional 
symmetric  key  ciphers,  and  the  recently  introduced  public-key  cryptomechanism 
characterized  by  the  ability  to  divulge  one  of  the  keys  in  the  matched 
(encryption,  decryption)  pair  without  leading  to  compromise  of  the  remaining 
key. 

The  development  hinges  on  two  mathematical  means  of  generalization, 
namely  (a)  the  transformation  of  the  message  space  from  Zp  to  Zp[x]/(g(x)) 

where  g(x)  is  any  irreducible  polynomial  over  Z  ,  the  finite  field  of 
integers  modulo  a  prime  number  p,  to  Zp[x]/(f(x})  where  f(x)  “IT  gjOO* 

the  product  of  distinct  Irreducible  polynomials  over  Zp;  and  (b^  the 
direct  summation  of  such  quotient  ring  terms,  isomorphically  realized  by 
invoking  the  Chinese  remainder  theorem  for  integers. 

The  newly  endowed  polynomial  structure  is  then  exploited  by 
concatenating  these  ciphers  with  a  data  scrambling  technique. 

Incorporation  of  the  resulting  cryptostructures  into  the  demanding 
environment  of  modem  multi-user  secure  communications  networks  is 
discussed  and  analyzed. 

Howard  M.  Shao,  Techniques  for  Signal  Processing  of  Image  Data,  Ph.D. 
Dissertation,  University  of  Southern  California,  May  1983. 

Abstract 

Techniques  for  processing  Images  have  become  of  increasing  interest  in 
recent  years.  Among  them,  very  often,  a  two-dimensional  filtering  operation 
has  been  required. 

The  commonly  used  FFT  method  to  compute  a  2-D  convolution  is  considered 
to  Involve  too  many  complex  multiplications  as  well  as  an  undesired  matrix 
transpose  operation.  New  algorithms  are  being  developed  to  cope  with  such 
problems . 
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It  was  first  shown  in  the  work  of  Nussbaumer  and  Quandalle  [l-3] 
that  a  certain  type  of  polynomial  transform  could  be  used  to  compute 
two-dimensional  convolutions.  Truong  and  Reed  [4]  modified  the  polynomial 
transform  with  the  Chinese  remainder  theorem  to  obtain  a  more  efficient 
algorithm.  By  using  a  new  factorization  of  the  modulus,  a  new  algorithm, 
called  the  Fast  Polynomial  Transform  (FPT)  was  developed.  This  technique 
takes  advantage  of  an  FFT-type  architectural  structure  to  compute  the  2-D 
convolution  more  rapidly  than  can  be  accomplished  with  the  standard  2-D 
FFT  structure. 

In  this  dissertation,  the  fast  polynomial  transform  is  further 
elaborated  and  implemented  on  a  general  purpose  computer.  It  is  shown  also 
how  to  realize  a  high  speed  array  processor.  Methods  to  cope  with  finite 
core  memory  limitations  are  prevented  and  actual  speed  performances  are 
obtained.  Since  50%  of  the  complex  multiplications  in  a  standard  2-D  FFT 
are  replaced  by  word  shifts  in  the  FPT,  experiments  show  that  the  FPT  has  a 
20%  speed  improvement  over  the  ordinary  2-D  FFT  for  computing  2-D  Convolutions. 

The  FPT  algorithm  developed  in  this  thesis  has  been  successfully 
applied  to  the  synthetic  aperture  radar  (SAR)  imaging  problem  [5] .  This 
involves  the  computation  of  two-dimensional  cross  correlation  of  radar  echo 
data  with  the  response  or  point  spread  function  of  a  point  target. 

Another  potentially  important  application  is  the  detection  of  a  moving 
target  in  optical  image  data  by  a  three-dimensional  space-time  filter  bank. 
Herein  an  optimal  3-D  filter  is  developed  for  detecting  optical  paths. 
Performances  is  computed  on  the  assumption  of  a  Markov  clutter  model.  The 
new  algorithm  shows  a  dramatic  increase  in  detectability  as  a  function  of 
target  velocity.  A  three-dimensional  FPT  algorithm  is  also  developed  to 
perform  the  3-D  space-time  filtering.  This  method  improves  the  computation 
speed  of  3-D  digital  filtering  significantly. 

Another  often  more  difficult  task  than  the  above  is  the  detection  of 
weak  stationary  targets  in  optical  image  data.  The  method  of  standard 
detection  is  often  unsatisfactory.  Thus  in  this  thesis  an  optimum  adaptive 
algorithm  is  developed  for  the  detection  of  a  signal  image  in  a  scene  with 
a  correlated  reference  scene  without  signal.  The  algorithm  is  implemented 
with  a  computer  program  utilizing  actual  image  data.  Experiments  show  that 
better  than  an  average  6  dB  lmprovment  can  be  achieved  over  the  method  of 
direct  correlation. 

For  real  world  applications  both  the  3-D  space-time  filtering  and  the 
optimum  adaptive  filtering  require  an  excessive  amount  of  computation.  The 
fast  polynomial  transform  and  other  similar  techniques  are  efficient 
computing  algorithms  for  improving  digital  signal  processing.  It  is  expected 
that  the  results  of  this  thesis  will  lead  to  the  development  of  more  practical 
and  useful  techniques  for  image  and  data  processing  problems. 


Chiunn-Shyong  Yeh,  Parallel  VLSI  Architectures  of  Signal  and  Communication 
Processors!  Ph.D.  Dissertation,  University  of  Southern  California,  June  1983. 

Abstract 

In  this  dissertation  VLSI  architectures  of  signal  and  communication 
processors  are  developed.  The  increasing  complexity  of  VLSI  devices  is 
a  challenge  to  both  architecture  and  circuit  designers.  The  development  costs 
and  turnaround  time  are  extremely  important  for  designing  custom  VLSI  devices. 
Herein  approaches  and  the  difficulties  associated  with  the  design  of  the  VLSI 
device  are  investigated.  The  properties  that  a  VLSI  architecture  should 
possess  in  order  to  minimize  the  design  costs  are  discussed.  The  important 
properties  of  VLSI  architectures  that  need  to  be  addressed  are  interconnects, 
regularity,  control,  modularity  and  concurrency. 

In  this  dissertation,  the  VLSI  architectures  of  three  signal  processors 
are  developed.  The  first  design  is  a  VLSI  architecture  for  the  arithmetic  of 
the  finite  field  GF(2m)  of  2m  elements.  Algorithms  are  derived  to  perform 
the  multiplication  and  inverse  in  the  finite  field  GF(2m) .  Then  two 
architectures  are  designed  to  realize  the  multiplication  algorithm.  These 
multipliers  then  are  used  to  design  an  inversion  circuit  for  the  finite  field 
GF(2m) . 

The  second  system  is  a  VLSI  architecture  for  digital  filters  which  uses 
the  Fermat  number  transform  (FNT) .  A  pipeline  structure  is  designed  to 
confute  a  128-point  FNT.  A  major  drawback  of  the  FNT  is  the  conflict  between 
long  transform  lengths  and  a  large  dynamic  range.  A  generalized  overlap-save 
algorithm  is  developed  to  solve  the  conflict  of  the  FNT.  Then  a  VLSI 
architecture  is  designed  to  realize  the  generalized  overlap-save  algorithm 
for  a  digital  filter  using  the  FNT. 

The  third  design  is  a  VLSI  architecture  for  digital  filtering  which  uses 
a  complex  number-theoretic  transform  (CNT) .  A  one-dimensional  array  is  designed 
to  compute  a  248-point  CNT .  Then  a  parallel  VLSI  architecture  is  designed  to 
realize  a  digital  filter  by  this  CNT  and  the  generalized  overlap-save  algorithm. 

The  architectures  of  these  signal  and  communication  processors  are  simple 
and  modular,  and  possess  the  properties  of  regularity  and  concurrency. 

Therefore  they  are  well-suited  for  implementation  in  VLSI  devices. 
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